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Abstract 

We present a study of B semileptonic decays into charm final states based on 211.7 fb~^ of data 
collected at the T(4S) resonance with the BABAR detector at the PEP-II e~^e~ storage ring. Using 
a novel technique based on the simultaneous fit of a set of variables reconstructed on the recoil 
of a B tagged in an hadronic decay mode, we measure the relative branching fractions r(B~ 
D^l-vi,)/T{B- DXe-Uf) = 0.210 ±0.017 (stat.) ±0.021 (syst.), r(B- ^ D*°i-i?i,)/T{B~ 
DXe-ue) = 0.611±0.022 (stat.)±0.027 (syst.) andr(B" ^ D**°{D^*K)£~ue)/T{B- DXrun) = 
0.173 ± 0.017 (stat.) ± 0.021 (syst.). 
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1 Introduction 



The determination of the individual exclusive branching fractions of B ^ X^lvi decays^ is im- 
portant for the study of the dynamics of semileptonic decays of the B meson. Precise data on 
the spectroscopy of the hadronic system is needed to be compared with theoretical predictions, 
obtained under different assumptions 012101 ID- This is especially important since the mass of 
the hadronic system, recoiling against the leptonic pair in the decay, is a crucial parameter in 
the extraction of \Vcb\ in exclusive semileptonic decays, in isolating B — > Xu^v^ from B — > Xcii^t 
decays for determining | Vub \ , and also in the extraction of the heavy quark masses and other non- 
perturbative Operator Product Expansion parameters from the distribution of spectral moments 
in semileptonic decays. This is exemplified by the fact that one of the leading sources of systematic 
uncertainty in the extraction of \Vcb\ from the exclusive decay B — > D*iv'£ is our limited knowledge 
of the background due to B ^ D*it£i'i 0. Reducing this uncertainty is of crucial importance, as 
lattice calculations promise to improve the theoretical accuracy on the form factor normalization 
to 2-3 %. Improvements in our knowledge of i? ^ Xcii'i decays will also benefit the accuracy in 
the extraction of \Vub\, as analyses are extending their probe into kinematical regions where these 
decays represent a sizable background. 

The first determination of the fractions of D** states in semileptonic b decays was obtained by 
the ALEPH |H] , DELPHI [71 |S] and OPAL experiments at LEP, where some information on the 
mass distributions of these states was also extracted, and by CLEO [1^ at CESR. In particular, the 
kinematics at LEP, where the B hadrons had a significant boost and were separated in opposite 
hemispheres, allowed analyses based on the topology of the particle tracks with respect to the 
decay vertices. On the other hand, primary fragmentation particles represented a background and 
statistics were limited. More recently, new results have been obtained by the DO [IV experiment 
at the Tevatron. 

Not only are the data samples at the B factories much larger than those obtained at earlier 
collider experiments, in addition the feasibility to fully reconstruct one of the two B mesons pro- 
duced exclusively in the decay of the T{4S) resonance permits the study of semileptonic decays of 
the other B meson in an almost unbiased way. This allows the study of exclusive semileptonic B 
decays involving charm mesons, specifically D, D*, resonant D** and non-resonant D^*h with the 
BABAR data. 

Measurements of the largest B meson branching fraction, B{B° D*-£+vi>), by the CLEO, 
the LEP and the BABAR and BELLE experiments, based on different reconstruction techniques, 
need to be reconciled. The current average value has a x^/ndof = 14.8/7 with a probability of 
just 3.8 %. 

Furthermore, there have been so far only very few studies of semileptonic decays to D** and 
non-resonant final states at the i?-factories. BELLE has recently reported the determination of the 
D^*^ttIvi decay branching fraction 13 . 

This study uses a novel technique to extract the exclusive branching fractions^ for B~ — > 
D^£~^i'£, B~ D*^£~i'£ and B~ D**^£~V£ {£ = e oi fi), where D**^ denotes here an hadronic 
final state containing a charm meson and with total mass above that of the D* state, therefore 
including both D* mesons and non-resonant IT states. 

^Here Xc refers to any charm state, Xu to any charmless particle. 

^Charge-conjugate modes are implied throughout this paper, unless explicitly stated otherwise. 
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2 The BaBAR Detector and Dataset 



This analysis is based on data collected with the BABAR detector at the PEP-II storage ring. 
The total integrated luminosity of the data set is 211.7 fb~^ collected on the T{4S). The corre- 
sponding number of produced BB pairs is 239 million. The BABAR detector is described in detail 
elsewhere [T^. Charged-particle trajectories are measured by a 5-layer double-sided silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH), both operating in a 1.5-T solenoidal magnetic 
field. Charged-particle identification is provided by the average energy loss (dE/dx) in the tracking 
devices and by an internally reflecting ring-imaging Cherenkov detector. Photons are detected by 
a CsI(Tl) electromagnetic calorimeter. Muons are identified by the instrumented magnetic-flux 
return (IFR). We use Monte Carlo simulations of the BABAR detector based on GEANT 15^ to vali- 
date the event reconstruction and the fitting technique used to extract the individual semileptonic 
branching fractions. 

3 Overview of the Analysis Method 

The technique developed for this analysis identifies a set of variables, which can be used to discrim- 
inate between the individual semileptonic decay modes. These variables are i) the missing mass 
squared reconstructed with respect to the D£ system, m^j^^, ii) the lepton momentum in the B 
rest frame and iii) the number of reconstructed tracks in addition to those used for reconstructing 
the D state and the lepton. In order to reduce the sensitivity to Monte Carlo modelling, shapes 
of the discriminating variables for the D^, D*^ and D**^ decays are extracted from data samples 
that are highly enriched with the individual decay modes. 

Semileptonic B decays containing one fully reconstructed D meson and recoiling against a fully 
reconstructed B in hadronic decay modes (Btag) are selected. We reconstruct Btag decays of the 
type B DY, where D refers to a charm meson, and Y represents a collection of hadrons with 
a total charge of ±1, composed of nivr^ + n2K^ + n^Kg + 71471*^, where ni + n2 < 6, < 3, and 
n4 < 3. Using and D*^ as seeds for B^ decays, we reconstruct about 1000 different decay 
chains. The kinematic consistency of a Btag candidate with a B meson decay is checked using 



two variables, the beam-energy substituted mass niES = \I^/^~Pb ™d the energy difference, 
/S.E = Eb — y/s/2. Here y/s refers to the total energy, and pB and Eb denote the momentum and 
energy of the Btag candidate in the T(45) center of mass frame. For correctly identified Btag, the 
mES distribution peaks at the B meson mass, while IS.E is consistent with zero. The signal region 
for tagging B candidates is defined as niES > 5.27 GeV/c^. 

The analysis exploits the presence of two charmed mesons in the final state: one used as a 
seed for the exclusive reconstruction of the Btag, and another in the B semileptonic decay, to 
obtain a high reconstruction efficiency. The reconstruction starts from the B semileptonic decay, 
identifying a charm meson of the correct charge-flavor correlation with a lepton with momentum 
in the T{AS) frame higher than 0.6 GeV/c and then selects the Btag with the smallest A.E among 
those which do not overlap with the reconstructed charm meson associated with the lepton. 
candidates are reconstructed in nine decay modes: — > ii'^vr"*", K^tt^tt^, K^-k^tt^-k' ^ K^tt^-k^ , 
K^tt'^ t:~ t:^ , K^tt^ , K~^K~ , tt^tt~ , K'^K^. candidates are reconstructed in seven decay modes: 
D+ K^7r+7r+, ivr-vr+TT+vrO, K^tt+ , K^tt+tt^, K+K-'k+, KIK+ , i^gvr+vr+vr-. 

This procedure is designed to maximize the efficiency of the Dil. reconstruction and, as a result, 
the Btag sample shows high purity, without need of further cuts. 

Then, B~ DXi^ug events, where X can be either nothing or any particle(s) from a charged 
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Figure 1: The niES distributions (data points with statistical errors) for the hadronic Bfag with 
B~ DXe~U(, (left) and B~ DXjjL~9^ (right) events. The solid line represents a fit to the total 
distribution, the dotted line the combinatorial background. The fitted signal yields are respectively 
2409 lb 122 and 1624 it 130 for the electron and muon samples. 



B semileptonic decay into an higher mass charm state (or a non resonant state), are identified by 
relatively loose selection criteria: we require the charm meson ground state invariant masses Mj^o 
and Mu+ to be in the range 1.85 < M^o < 1.88 GeV/c^ and 1.853 < < 1.883 GeV/c^ and 

the cosine of the angle between the D candidate and the lepton in the T{AS) frame to be less than 
zero, to reduce background from non B semileptonic decays. After these loose selection criteria, 
the sample contains leptons from prompt and cascade B meson decays (i.e. the lepton does not 
come directly from the B) plus various background sources (e.g. photon conversions and tt*^ Dalitz 
decays, combinatoric BB background) that need to be subtracted. The contamination from cascade 
B meson decays (about 15% of the total sample) is removed by using the simulated Monte Carlo 
distributions for these backgrounds. These events are reweighted to account for possible differences 
between the branching fractions used in our Monte Carlo simulation and the latest experimental 
measurements. The photon conversion and vr'^ Dalitz decay background (less than 2% of the total 
electron sample) is removed by using a dedicated algorithm, which performs the reconstruction 
of vertices between tracks of opposite charge whose invariant mass is compatible with a photon 
conversion or a vr'^ Dalitz decay. 

Figure n shows the tues distribution for the Btag candidates in the B^ DXi^ui sample. 
We fit this distribution with empirical functions: a Crystal Ball ^01 for the signal and an AR- 
GUS function |17j for the background. Combinatoric BB background is removed by performing 
a sideband subtraction by scaling the number of background events in the Btag fnES sideband 
region 5.21 < mEs < 5.26 GeV/c^ to the integral of the background function in the signal region. 
Cross- feed effects {B^^g candidates erroneously reconstructed as a neutral B) are corrected for using 
Monte Carlo simulation. We estimate the fraction of B^ events in the reconstructed B~ sample to 
be 7%. 

To determine the total reconstruction efficiencies for individual signal decays using Monte Carlo 
simulation, the fitted yield of events in the niES distribution of the Btag decays is compared to the 
produced number of semileptonic decays in tagged events. 

The inclusive distributions of the discriminating variables are determined for the selected B~ 
DXi~v'i samples. The missing mass squared in the event is defined as: 
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mLss = (PiT) - p{Btag) - p{D^^+) - p{fi)f (1) 

in terms of the various particle four-momenta and for correctly measured signal events it corresponds 
to the invariant mass of the Xv^ system in the B~ D^'^Xl^'ui decay. Thus for B~ D^£~i'i 
decays, the missing mass squared corresponds to the neutrino mass and peaks at zero. For semilep- 
tonic decays to excited D mesons or non-resonant states, the distribution is shifted towards higher 
masses. The selection of fully reconstructed events results in an excellent missing mass resolution 
of 0.04 GeV^/c^, an order of magnitude lower compared to non-tagged analyses [TS] . 

The lepton momentum in the T(4S') frame is the second variable which can be used to some 
extent to discriminate among the different B semileptonic decays. Leptons directly produced from 
B decays will have momenta greater than secondary cascade decays of the type B ^ D ^ i. For 
true B semileptonic decays, the lepton momentum is peaked between 1.5 and 2 GeV/c. Leptons 
from B Dlui will usually have lower momenta than B D*£i'i, because for pD — > 0, the decay 
B Diu£ is suppressed, so D mesons will have higher momenta. Leptons from B D**{D^*^7r)ii'i 
will have much lower momenta, due to the reduced phase space available. 

The multiplicity of reconstructed charged tracks in the event, in excess of those accounted for 
the Btag recoil side, the lepton and the D^'~^, offers as the missing mass a very high discriminating 
power to separate the different B semileptonic decays. Every charged track is required to have its 
point of closest approach to the interaction point less than 10 cm along the beam axis and less 
than 1.5 cm transverse to the beam axis. The number of additional tracks should be zero for pure 
B — > D^'~^£i'£ decays and larger than zero for those decays involving heavier charm states, provided 
these decay into a lighter D accompanied by charged particles. 

Figure 121 shows the reconstructed missing mass squared, lepton momentum and additional 
charged track multiplicity distributions for a Monte Carlo sample of B^ DXi~ui events. The 
different shapes of the B~ D^*'**^^i~i?£ components are clearly visible. As it can be observed, a 
fraction of B~ D^*^'^l~i>i events shows a multiplicity of additional charged tracks greater than 
zero, coming from badly reconstructed tracks not used for the Btag reconstruction. 

4 Reconstruction of the Exclusive Decays B D^* z/£ 

To extract the relative branching ratios of the B^ D^*'**^^£^i>£ decays from the inclusive B~ — 
DXi~v'i sample. Probability Density Functions (PDFs) corresponding to the different components 
present in the inclusive distributions must be built. These PDFs are determined using data sets 
enriched in the different exclusive B^ D^*'**^^£~U£ decays and then parameterized as analytical 
functions, with the exception of the charged particle track multiplicities which, owing to their 
discrete nature, are best described by histograms. 

The selection of these enriched samples is optimized in terms of their purity, to introduce 
minimum biases with respect to the true underlying distributions. Background contributions in 
the exclusive samples include other B semileptonic decays, where one particle is not reconstructed 
or erroneously added to the charm meson candidate (feed-down and feed- up). We use a sample of 
B~ — > D^'~^X£~i'£ or B~ — > D*^'~^Xi~Ui events to reconstruct the exclusive decays. and 
candidates are selected as above. D* mesons are reconstructed by combining a D candidate and 
a pion in the D*~^ D^ir^(D^TT^) and D*^ D^tt^ decays. These D* candidates are required 
to have a mass difference Am = m/^^ — m/j within ±2,3,4 MeV/c^ of its nominal values for the 
decays D*~^ — > Z^+tt^, D*^ — > D^-k^ and D*^ —> D^Tr~^ respectively. For the decays with a photon. 
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Missing IMass Squared (GeV^/c") Additional Tracl< IMultiplicity 
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Lepton lUomentum (GeV/c) 



Figure 2: Distributions of the missing mass squared, the multiphcity of additional charged tracks 
and the lepton momentum spectrum for a Monte Carlo sample of B~ DXi^u^ events used for 
the relative branching ratio V{B~ — > D^'*^'**^l~U()/T{B~ — > DXl~v^) measurement; the different 
components are superimposed with different colors. 

D*^ D^'y, we require the mass difference Am = m£)j — rriD within ±9 MeV/c^ of its nominal 
value. D**'^ states or non resonant B~ D^*^iTi^v'i decays are reconstructed by combining a D or 
D* candidate with a charged pion. Neutral pions are not used at this stage of the reconstruction 
to avoid large combinatoric background. 

Using the four-momentum of the reconstructed candidate, we compute the corresponding 

missing mass squared . „r* in the event, defined as: 

<^ss,D(*•*'^ = (piT) - p{Btag) - p{D<**)) - p{i)f . (2) 

For true B^ — > D^*'**^^£~i'i semileptonic decays, this missing mass squared will peak at zero. 
Exclusive samples are then selected by applying the relevant cuts on this variable to select or 
remove, if we want to eliminate feed-down components, a specific B semileptonic decay. 

To select D*^i~i>£ decays, the variable "T-^iss D*o required to be in the range —0.35 < 
''^miss D*o ^ ^-^ GeV^/c^, while the unwanted decays are removed, by vetoing the peaking re- 
gion in the corresponding missing mass squared distribution, i.e. f^l^iss D**+ ^ ^'^ GeV^/c^ or 
<.^ss,D»^ < -0-6 GeVVc^ and m^^^^ ,,..0 < -0.2 GeVVc^- 

The exclusive decay B~ — > D^i~ui is selected by removing feed-down from D* and D** events. 
We require Tn^^gg £)»o < —0.05 GeV^/ to reduce feed-down from D*'^ — > D^y decays, n^'^^gg £)*o < 
-0.4 GeVVc^ or . „.o > 1-0 GeVVc"^ for D*° L»°7r° decays. Feed-down from D*+ events 
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is reduced by requiring I'^n'^iss D*+\ ^ ^■'^ GeV^/c^. Feed-down from D** events is reduced by 
requiring |m^.^^^^..o| > 0.8 GeVVc^ or m^.^^^^„o < -0.5 GeVVc^ for D**^ D*+tt- and 
£>**° ^ D*°7r° decays and ml^.^^ j^,,+ < -0.5 GeVVc^ for D**+ events. 

Exclusive B~ — > D**^i^V£ events are selected by requiring —0.25 < i^^miss D**+ ^ 0.85(0.25) 
GeVVc^ for D**° D+{D*+)tt~ decays. 

Non B semileptonic decay background and combinatoric BB background are subtracted as in 
the inclusive distribution. Table Q lists the composition of the reconstructed exclusive samples in 
terms of purity, feed-up (-down) fractions and non B semileptonic background. 



Table 1: Purity and feed-down(up) contributions for the reconstructed exclusive decay modes. 



Decay Mode 


Purity 


Feed-up (D) 


Feed-up (D*) 


Feed-down (D*) 


Feed-down (D**) 


background 




75% 






19.3% 


3.6% 


2.5% 


B- D*'^e-Pe 


91% 


3.7% 






2.8% 


2.6% 


B- D**"{D*+n')e-Ui 


85% 


0.27% 


5.5% 






9.4% 


B- D'""{D+n-)e-ui 


85% 


1.2% 


4.2% 






9.7% 



This reconstruction strategy has been validated to assess whether the exclusive selection proce- 
dure induces biases on the shapes of the discriminating variables. This has been done by comparing 
the distribution of these variables using different selection criteria on Monte Carlo simulated events 
of the different decay modes B~ D^*'**^^£~i>i. The agreement of these distributions is good: 
residual discrepancies are evaluated as systematic effects. 

The final step is the determination of the fractions of D^, D*^ and D**^ decays in the selected 
inclusive semileptonic sample, by a simultaneous fit to the missing mass squared, additional 
charged track multiplicity and lepton momentum distributions where the relative fractions are 
treated as free parameters together with the parameters describing the shapes of the discriminating 
variables. The PDFs for the three different fit components are extracted by fitting simultaneously 
the corresponding exclusive and inclusive reconstructed distributions. The missing mass squared 
PDF is the sum of three Gaussians, with independent means and widths, to account for the upper 
tail, which is mainly due to feed-down from B~ — > D*^£~r'i decays where the soft pion or the 
photon from the D*^ is missing. The D*^ [D**^) missing mass squared PDF is built by the product 
of a Gaussian (an exponential) with a polynomial function. The three D^, D*^ and D**^ lepton 
momentum spectra are parameterized by the product of a polynomial and an exponential which 
describes the low momentum tail of the spectra. Because we do not reconstruct D**^ states involving 
neutral pions, we use the B^ D*^l~i>i additional charged track multiplicity distribution to model 
B~ — > D**^£~r'i events, where the D**^ decays in neutral pions. 

The total PDF used to model the inclusive distribution is obtained by the sum of the three 
exclusive PDFs, whose relative weights correspond to the relative fractions of D^, D*^ and D**^ 
decays in the sample. The sum of the relative weights is constrained to be equal to one. This results 
in a 33-parameter fit which ensures that statistical correlations are properly taken into account and 
the uncertainties on the exclusive shapes, obtained on much reduced statistics compared to the 
inclusive sample, are correctly propagated into the statistical uncertainties on the D^, D*^ and 
D**^ fractions. The fit also takes into account the feed-up and feed-down fractions in the exclusive 
shapes, obtained from simulation and contributing to the systematic uncertainty. 

The fit results for the B~ ^ DXl^vi^ are shown in Figure IHl and IH Figure 01 shows the fit 
result for the exclusive missing mass squared and lepton momentum distributions, which are used 
to build the PDF used to fit the corresponding inclusive distributions, shown in Figure IH The 
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simultaneous fit has a total x^/ndof = 237/205. 



5 Systematic Studies 

Different sources of systematic uncertainties have been estimated. These are divided in various 
categories. Detector related systematics may arise from discrepancies in the simulation modelling 
of the track reconstruction and efficiency, particle identification and neutral particle reconstruction. 

Uncertainties related to the reconstruction of charged tracks are determined by evaluating the 
fit stability using different track selection criteria and by a toy Monte Carlo study in which we 
fluctuate the track multiplicity bin content according to the tracking efficiency uncertainty. The 
systematic error due to the reconstruction of neutral particles is studied by varying the resolution 
and efficiency to match those found in control samples on data. We estimate the systematic error 
due to particle identification by varying the electron and muon identification efhciencies by 2% 
and 3%, respectively. The misidentification probabilities are varied by 15% for both electrons and 
muons. 

Different systematic uncertainties are evaluated for the inclusive and exclusive reconstruction 

procedure. For the inclusive reconstruction of the B~ DXI'~ui samples, the weighting factors 
used to subtract B cascade decays are varied within their errors and the differences in the fitted 
relative branching ratios are included in the systematic uncertainty. We also evaluate systematic 
effects associated to the misidentification rate and the tracking efficiency of the photon conversion 
and vr'' Dalitz decay reconstruction algorithm. The uncertainty of the Btag combinatorial back- 
ground subtraction is estimated by evaluating differences in the shapes of this background in the 
sideband and in the signal region using Monte Carlo predictions. We also evaluate the dependence 
of the sideband scaling factors on the lepton momentum. The systematic error due to the un- 
certainty in the amount of flavor cross-feed is computed by varying its fraction by a conservative 
30%. The effect of uncertainty on reconstruction efficiency mostly cancels in the relative branching 
ratios; the uncertainty due to the finite Monte Carlo statistics and the niES fit we use to compute 
the relative reconstruction efficiency is included accordingly in the systematic error evaluation. 

For the exclusive reconstruction of B~ D^*'**')^i~iy£ decays, systematic uncertainties arise 
from the evaluation of the purity of each sample and residual biases on the reconstructed dis- 
tributions due to the selection criteria. Possible mismatching between data and Monte Carlo is 
accounted for by varying the feed-down and feed-up fractions within their errors and by evaluating 
as systematic uncertainties differences in the fit results. We also evaluate possible biases due to the 
use of reconstructed exclusive distributions to build the PDFs by replacing them with Monte Carlo 
signal distributions and performing the fit using these shapes. A separate systematic uncertainty 
is computed for the additional track multiplicity distribution by varying the relative ratio of 
D**^ D^*^^Tr~ and D**^ D^*^^Tr'^ events which is predicted by our Monte Carlo simulation. 
The uncertainty of the Btag combinatorial background subtraction and the cross-feed correction is 
computed as for the inclusive reconstruction. 

The fit performance has been extensively tested using simulated samples with varying fractions 
of the different decay modes. These tests show that the procedure adopted in this analysis is able 
to extract the decay fractions without significant biases; the statistical errors obtained from the 
fit are corrected according to the results on independent toy Monte Carlo tests and increased by 
15(20)% for the d'^*^^{D**^) component. We evaluated the systematics due to the fit technique 
by using different parameterizations of the exclusive shapes. Additional cross-checks were made, 
repeating the simultaneous fit of the exclusive shape parameters and the exclusive yields in different 
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Figure 3: Missing mass squared and lepton momentum distributions for the three selected exclusive 
decay modes, B~ D^£~i'£, B~ D*^i~r'£ and B~ D**^l~9^. The data are compared to 
the result of the overall fit (solid line) to the distributions for the selected exclusive and inclusive 
samples. 
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Missing lUlass Squared (GeV^/c") Additionai Number of Cliarged Traci<s 




Lepton Momentum (GeV/c) 

Figure 4: Distributions of the missing mass squared, the multiphcity of additional charged tracks 
and the lepton momentum spectrum for the selected B~ — > DXi'v^ events. The data are compared 
to the results of the overall fit (solid line). The PDFs for the three _d(*'**)o components are shown 
in different colors, from lighter [D**^) to darker {D*^). 
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configurations, e.g. by fixing the PDF parameters for the exclusive distributions or by only fitting 
the inclusive distributions. 

6 Results 

The extraction of the relative fractions of D^, D*^ and D**^ decays in the selected inclusive semilep- 
tonic sample gives the relative branching ratios T{B- D^*'**^°e-Ui)/T{B- DX£-ui). 

The measured branching ratios are reported in TablelSJand the different systematic uncertainties 
are shown in Table El 

Table 2: Fitted branching ratios T{B- D'^'*'^'**°i'i?e)/T{B- DXl'v^) with the corresponding 
statistical uncertainty. 



Branching Ratio 


Fit Results 


V{B- ^D"l-ui) 

r(B-^Dxe-Pf) 


(21.0 ± 1.44) % 


r{B-^D*"e-Pf) 

r(B-^DXf~ut) 


(61.1 ± 1.94)% 


r{B-^D**"e-pe) 

V{B-^DXl-ut) 


(17.3 ± 1.44) % 



Table 3: Systematic uncertainties in the measurement of r(i? ^ i'i)/T{B DXi vg). 





Systematic uncertainty on T(B D 


'■•"'>oe-pg)/r{B- ^Dxe-pi) 




B~ D^lvi 


B- D*^£ue 


B- D**^£ue 


Tracking efficiency 


0.009 


0.008 


0.004 


Neutral reconstruction 


0.001 


0.003 


0.0009 


Electron ID 


0.0008 


0.002 


0.0007 


Muon ID 


0.006 


0.02 


0.005 


Inclusive Reconstruction 




Cascade decay background 


0.01 


0.01 


0.01 


Conversion and Dalitz decay background 


0.001 


0.004 


0.001 


MC statistics and Mes fit 


0.012 


0.01 


0.011 


Cross-feed corrections 


0.001 


0.002 


0.004 


Exclusive Reconstruction 




Feed-down and feed-up corrections 


0.003 


0.002 


0.002 


MC signal shapes 


0.002 


0.002 


0.007 


D**^ track multiplicity 


0.005 


0.002 


0.005 


Mes fit 


0.0007 


0.004 


0.003 


Cross-feed corrections 


0.0003 


0.0006 


0.0009 


Fit Technique 




Missing mass and lepton momentum PDF 


0.007 


0.002 


0.008 


Total Systematic Error 


0.021 


0.027 


0.021 



No direct measurement of the relative branching ratios T(B^ D^*'**^'^£^ i^g) /T{B^ DXl^v^) 
is present in the literature. In the absence of any measurement [5j for B semileptonic decays into 
charm final states Xc other than D meson states, we can assume the relation: 
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r{B- XJ-ue) r{B- DXl-vi) ^ ' 

which can be used with the current best values ^1] for Xli'i, b ulvi and the branching 

fractions B~ D^*'**')^£i>£ to make a comparison with our results. 

The -B^ — > D**^l^i'£ component is also sensitive to the presence of B^ D^*^mTi^ v'l events 
in the inclusive sample B~ — > DX despite we do not explicitly introduce in the fit the 
corresponding PDF. By floating its parameters in the fit, the D** PDF can model any excess of 
events in the inclusive sample B~ — > DXi~vi, which can be associated to a small B~ D^*^mT£~i?i 
component. 

7 Conclusions 

Preliminary results of the relative branching fractions of individual semileptonic B decays have 
been obtained for the B~ decaying into D^£~U£, D*^£~U£ and D**^l~V£. These measurements are 
based on a sample of 239 million T(45') BB decays collected by the BABAR experiment at the 
PEP-II asymmetric-energy B factory at SLAC, in events in which one B meson decaying to a 
hadronic final state is fully reconstructed. A novel technique, based on a multi-parametric fit to 
a set of discriminating variables in an inclusive sample of B~ — > DXi~U£ events, has been used 
to measure the relative branching fractions for B~ D^,D*^ and D**^{D^*'^7r)£~i>£ decays. In 
order to reduce the sensitivity to Monte Carlo modelling, which is especially limited by the poor 
knowledge of the B~ — > D**^ {D^*^ ■7r)£~ i>£ decays, the shapes of the discriminating variables for the 
D^, D*^ and D**^ decays are extracted from selected data samples of exclusive decays. 
Correcting the statistical errors as discussed in section 5, we obtain: 

0.210 ± 0.017 (stat.) ± 0.021 (syst.) 
0.611 ± 0.022 (stat.) ± 0.027 (syst.) 
0.173 ± 0.017 (stat.) ± 0.021 (syst.) 

The BELLE result jl3j is consistent with our measurement. The extension of this technique 
to the B^ decay modes wih increase the accuracy of the branching fraction mea- 

surements of B semileptonic decay in D meson final states and will help the understanding of the 
observed difference between the rate for the inclusive semileptonic decay and the sum of the rates 
for the inclusive modes [Sj. Moreover, an accurate measurement of the B D**£v branching 
fraction, which is made possible by the use of this method, is important to understand the decay 
mechanism of this channel, which provides crucial tests of the Heavy Quark Effective Theory and 
the QCD Sum rules [El. 
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